Onset of classical behaviour after a phase transition 
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We analyze the onset of classical behaviour in a scalar field after a continuous phase transition, 
in which the system-field, the long wavelength order parameter of the model, interacts with an 
environment, of its own short-wavelength modes. We compute the decoherence time for the system- 
field modes from the master equation and compare it with the other time scales of the model. 
Within our approximations the decoherence time is in general the smallest dynamical time scale. 
Demanding diagonalisation of the decoherence functional produces identical results. The inclusion 
of other environmental fields makes diagonalisation occur even earlier. 
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I. INTRODUCTION 

The standard big bang cosmological model of the early 
universe assumes a period of rapid cooling, giving a 
strong likelihood of phase transitions, at the grand uni- 
fied and electroweak scales in particular. What inter- 
ests us in this talk is the way in which phase transitions 
naturally take us from a quantum to classical description 
of the universe. 

That (continuous) transitions should move us rapidly 
to classical behaviour is not surprising. Classical be- 
haviour has two attributes: 

• Classical correlations: By this is meant that the 
Wigner function(al) VF[7r,(/)] peaks on classical 
phase-space trajectories, with a probabilistic inter- 
pretation. 

• Diagonalisation: By this is meant that the density 
matrix p{t) should become (approximately) diag- 
onal, in this case in a field basis. Alternatively, 
we can demand diagonalisation of the decoherence 
functional. In either case a probabilistic description 
(no quantum interference) is obtained. 

From the papers of Guth and Pi j^j] onwards, it has 
been appreciated that unstable modes lead to correla- 
tions through squeezing. On the other hand, we under- 
stand diagonalisation to be an almost inevitable conse- 
quence of tracing over the 'environment' of the 'system' 
modes. 

Continuous transitions supply both ingredients. 
Firstly, the field ordering after such a transition is due 
to the growth in amplitude of unstable long-wavelength 



modes, which arise automatically from unstable maxima 
in the potential. Secondly, the stable short-wavelength 
modes of the field, together with all the other fields with 
which it interacts, form an environment whose coarse- 
graining enforces diagonalisation and makes the long- 
wavelength modes decohere. 

While there are few doubts about the classical out- 
come, to quantify these general observations is difficult 
because, with fields, we are dealing with infinite degree 
of freedom systems. One of us (F.L) has shown elsewhere 
3] how classical correlations arise in quantum mechani- 
cal systems that mimic the field theory that we shall con- 
sider here, and we refer the reader to that paper for the 
role that classical correlations play. Our concern in this 
talk is, rather, with diagonalisation, equally necessary for 
the onset of classical behaviour. This is determined both 
through the master equation for the evolution of the den- 
sity matrix and the decoherence functional, whose role is 
to describe consistent histories. 

This talk builds upon earlier work by us 13, Isl, ld| , and 
we refer the reader to it for much of the technical details. 
We restrict ourselves to flat space-time. The extension 
to non-trivial metrics is straightforward in principle. See 
the talk of Lombardo in these same proceedings , which 
complements this. 

II. BASIC IDEAS 

The evolution of a quantum field as it falls out of equi- 
librium at a transition is determined in large part by its 
behaviour at early times, before interactions have time 
to take effect. To be concrete, consider a real scalar 
field described by a Z2-symmetry breaking action 

(m' > 0) 
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with symmetry breaking scale 77^ — 6/i^/A. On heating, 
this shows a continuous transition, with critical temper- 
ature = 2?7^. If, by virtue of the expansion of the uni- 
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verse the system is very rapidly cooled (quenched) from 
T > Tc to r < Tc, the initial stages of the transition can 
be described by a free field theory with inverted mass 
—H^ < 0. The state of the field is initially concentrated 
on the local maximum of the potential, and spreads out 
with time. This description is valid for short times, until 
the field wave functional explores the ground states of 
the potential. 

The (/)-field ordering after the transition is due to 
the growth in amplitude of its unstable long-wavelength 
modes, which we term <f>^(x). For an instantaneous 
quench these have wave-number k < fj, for all time. 
Although the situation is more complicated for slower 
quenches, until the transition is complete there are al- 
ways unstable modes. As a complement to these, we 
anticipate that the stable short-wavelength modes of the 
field (/>>(a;), where 

0(x) = (t)<{x) + 4>>{x), 

will form an environment whose coarse-graining makes 
the long- wavelength modes decohere In practice, the 
boundary between stable and unstable is not crucially 
important, provided there is time enough for the power in 
the field fluctuations to be firmly in the long-wavelength 
modes. This requires weak coupling A ^ 1. Of course, 
all the other fields with which (f) interacts will contribute 
to its decoherence, but for the moment we ignore such 
fields (before returning to them in the last section). 
After splitting, the action Q can be written as 

=5[0<]+5[(/.>]+5int[</'<,<A>], (2) 

where the interaction term is dominated 0, |^ by its bi- 
quadratic term 

1 



S'int [4>< 
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A / d^x (j)'^{x)<j)iix). 



(3) 



The total density matrix (for the system and bath 
fields) is defined by 

and we assume that, initially, the thermal system and its 
environment are not correlated. 

On tracing out the short-wavelength modes, the re- 
duced density matrix 

whose diagonalisation determines the onset of classical 
behaviour, evolves as 



= J d(j}^i J dcj)^^ Jr[t,tt] Pr[tt], 

where J^[t,ti] is the evolution operator 



Mt,t,] 



Vc^< P0< exp{i5cG[0^,0<]}. 



S'cg[0o '?^<] is the coarse-grained effective action, of the 
closed time-path form 

All the information about the effect of the environment 
is encoded in ^S'[(/)^, 0^] through the influence functional 
(or Feynman- Vernon functional 0) 

F[0+,(/)<]=exp{z<55[0+,0<]}, 

giving 6S a well defined diagrammatic expansion. 



III. THE MASTER EQUATION 

To see how the diagonalisation of occurs, we con- 
struct the master equation^ which casts its evolution in 
differential form. As a first approximation, we make a 
saddle-point approximation for J,, in (0J , 

J^[(l4 ,tf\(f)f ,(1)7 ,ti] Ki e-x:p{iScG[(l)tv(l)c\\)^ (5) 
In (|3J| (fy^-y is the solution to the equation of motion 



SReScG 



= 0, 



(4) 



with boundary conditions 0^(io) = a-nd ip^^it) = (f>f . 

It is very difficult to solve this equation analytically. 
We exploit the fact that, after the transition, the field 
cannot be homogeneous in one of its groundstates 4> = rj 
or (f) — ~rj because of causality 10]. As a result there 
is an effective 'domain' structure in which the domain 
boundaries are 'walls' across which (j) ffips from one 
groundstate to the other. Further, these domains have a 
characteristic size ^, where = Trfco labels the domi- 
nant momentum in the power of the (/)-field fluctuations 
as the unstable long-wavelength modes grow exponen- 
tially. For simplicity, we adopt a 'minisuperspace' ap- 
proximation, in which we assume regular domains, en- 
abling (pcxix, s) to be written as 

0ci(x,s) = /(s,t)<i>(x)$(y)$(z), (6) 

where $(0) = <I>(^) = 0, and 

$(a; + C) = -*(x). 

f{s,t) satisfies f{0,t) = 0i and f{t,t) = (pf. We write it 
as 

f{s,t) = (j)iUi{s,t)+(j>{U2{s,t). (7) 

In we made the simplest choice for ^{x), 

^{x) — sin fco^;. 

Extensions to include more Fourier modes are straight- 
forward in principle, but our work in |^ was sufficient to 
show that the results only depend weakly on the details 
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of the domain function for few Fourier modes. In 

the Hght of the more qualitative comments made here, 
we refer the reader again to @ for details. On the other 
hand, the Ui{s,t) are solutions of the mode equation 
for wavenumber fco during the quench, with boundary 
conditions ui(0,t) = 1, ui{t,t) — and U2(0, i) = 0, 

U2{t,t) = 1. 

In order to obtain the master equation we must com- 
pute the final time derivative of the propagator J^. After 
that, all the dependence on the initial field configura- 
tions (jyf^ (coming from the classical solutions 0^) must 
be eliminated. Assuming that the unstable growth has 
implemented diagonalisation before back-reaction is im- 
portant, Jr can be determined, approximately, from the 
free propagators as 



strongly on the properties of the environment. It satisfies 



to 



ds D{ko, s), 



(10) 



and corresponds to the time after which we are able to 
distinguish between two different field amplitudes, inside 
a given volume V. 

To terms up to order and one loop in the h expan- 
sion (we continue to work in units in which fi, = fc^ = 1), 
the influence action due to the biquadratic interaction 
between system and environment has imaginary part 



ImJS* = 



d^x / d%A{x)N{x,y)A{y), (11) 



Jo[t, ti 



''l?0< / "'l?0<exp{i[5o(0+)-S'o(r)]} 



(8) 

where 5*0 is the free-field action. This satisfies the general 
identities B 



A±r ( +^ "2(i,i) , ui[s,t) 



ui{t,t) 



Ul{t,t) 



Jo 



which allow us to remove the initial field configurations 



^ and obtain the master equation. 



Even with these simplifications the full equation is very 
complicated, but it is sufficient to calculate the correc- 
tion to the usual unitary evolution coming from the noise 
(diffusion) kernels (to be defined later). The result reads 



ip, = (0+f|[i/,Pr]|0<f) - tVA^D{ujo,t)p, + ... (9) 
where D is the diffusion coefficient and 



</>r2)/2 



for the final field configurations (henceforth we drop the 
suffix). The ellipsis denotes other terms coming from the 
time derivative that do not contribute to the diffusive 
effects. V is understood as the minimal volume inside 
which there are no coherent superpositions of macroscop- 
ically distinguishable states for the field. 

The effect of the diffusion coefficient on the decoher- 
ence process can be seen by considering the following 
approximate solution to the master equation: 



exp 



FA^ / ds D{ko,s) 



where is the solution of the unitary part of the master 
equation (i.e. without environment). The system will 
decohere when the non-diagonal elements of the reduced 
density matrix are much smaller than the diagonal ones. 

The decoherence time to sets the scale after which we 
have a classical system-field configuration, and depends 



where N{x,y) = jX^RcG^^{x,y) is the noise (diffusion) 
kernel and G^_^_{x,y) is the thermal short- wavelength 
closed time-path correlator. 

Explicit calculation shows that D{kQ, t) takes the form 



where 



m(s, t) 



'o,t)= I 
Jo 



ds w(s, t) F{ko, s, t) 



(12) 



U2{s,t) 



Ui{t,t) 



ui(s,t) 



and F{kQ, s, t) is built from the spatial Fourier transforms 
of the overlap of the diffusion kernel with the field profiles 
$(a;)$(y)$(z). 

In the integrand of H12|l u{s,t) is rapidly varying, 
driven by the unstable modes, and F{ko,s,t) is slowly 
varying. For long-wavelengths fco <C /i we have, approxi- 
mately, 



whereby 



F{ko,s,t)^OiN{ko^O;t-s)), 



D(fco,t) « F(fco,0,i) / dsu{s,t). 

Jo 



(13) 



That is, the diffusion coefficient factorises into the 
environmental term F, relatively insensitive to both 
wavenumber and time, and the rapidly growing integral 
that measures the classical growth of the unstable system 
modes that are ordered in the transition. 

To be specific, we restrict ourselves to the simplest 
case of an instantaneous quench from a temperature T = 
0{Tc) > Tc, for which 

smh[uJo{t - s)] sinh(^os) 

'^i = ^TT — 7\ — ' ""2 = ^-TT — 7^, (14) 

smh(cL'oi) smh[uJot) , 

where uS^ — p? — k^ p^ . It follows that 

u(s, i) = cosh^[a;o(i — s)], (15) 
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from whose end-point behaviour at s 
we find the even simpler result 



of the integral 



D{ko,t) ^ fi-^F{ko,0,t) u{0,t) 



{XTjA^nf exp[2^t], 
(16) 

assuming jito ^ 1- 

For more general quenches growth is more complicated 
than simple exponential behaviour but a similar separa- 
tion into fast and slow components applies. 

We have omitted a large amount of complicated techni- 
cal detail (see @): to give such a simple final result. This 
suggests that we could have reached the same conclusion 
more directly. 

We now indicate how we can obtain the same results 
by demanding consistent histories of the 4> field. 



IV. THE DECOHERENCE FUNCTIONAL 

The notion of consistent histories provides a parallel 
approach to classicality. Quantum evolution can be con- 
sidered as a coherent superposition of fine-grained histo- 
ries. If one defines the c-number field (/'(x) as specifying 
a fine-grained history, the quantum amplitude for that 



(we continue to work in units in 



history is ~ e*'^'" 
which h = 1). 

In the quantum open system approach that we have 
adopted here, we are concerned with coarse-grained his- 
tories 



(17) 



where a[(/)] is the filter function that defines the coarse- 
graining. 

From this we define the decoherence function for two 
coarse-grained histories as 



(18) 

a ] does not factorise because the histories (fr^ are 
not independent; they must assume identical values on a 
spacelike surface in the far future. 

A necessary and sufficient condition for the validity 
of the sum rules of probability theory (i.e. no quantum 
interference terms) is |Tll | 



ReV [a+,( 



0, 



(19) 



when ^ (although in most cases the stronger 
condition I?[a+,a^] « holds 12J). Such histories are 
consistent 13] . 

For our particular application, we wish to consider as 
a single coarse-grained history all those fine-grained ones 
where the full field remains close to a prescribed clas- 
sical field configuration (pci- The filter function takes the 
form 



[</)] ^ VJe 



Ctcl 



[J]. 



(20) 



In principle, we can examine general classical solutions 
for their consistency but, in practice, it is simplest to 
restrict ourselves to solutions of the form ©. In that 
case, we have made a de facto separation into long and 
short-wavelength modes whereby, in a saddle-point ap- 
proximation over J, 



As a result, 



\V{ 



cl ' ^c\l 



exp{-Im^5[0+ 0-]} 



(21) 



(22) 



For the instantaneous quench of Q, using the late time 
behaviour 0^ ^ e^^cj)^, Im5S'[(/)^, takes the form 

T/A2 /■* /■* 

IuySS — ds ds'e^''' e^f"' F{ko,s,s'). (23) 

fJ- Jo Jo 

From this viewpoint adjacent histories become consis- 
tent at the time tu, for which 



dt Im SS. (24) 



V. THE DECOHERENCE TIME 

We have used the same terminology for the time tjj 
since, on inspection, (|24|l is identical to H10() in defining 
the onset of classical behaviour. As we noted, in practice 
the use of the decoherence functional looks to be less re- 
strictive than the master equation, and we hope to show 
this elsewhere. 

For the moment what is of interest is whether to , based 
on linearisation of the model, occurs before backreaction 
sets in, to invalidate this assumption. When all the de- 
tails are taken onto account, whether from ^ or 114|l . 
tu satisfies 



exp(4^iD) 



or, equivalently 

exp(4/itD) = O 



(25) 



(26) 



For the rapid quenches considered here, linearisation 
manifestly breaks down by the time t* , for which {(l>'^)f ~ 
rj^, given by 



exp(2Mt*)=0f^ 



(27) 



The exponential factor, as always, arises from the growth 
of the unstable long-wavelength modes. The factor 
comes from the coth(/9a;/2) factor that encodes the initial 
Boltzmann distribution at temperature T >Tc. 

Our conservative choice is that the volume factor V 
is 0{fi~^) since fi~^ (the Compton wavelength) is the 
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smallest scale at which we need to look. With this choice 
it follows that 

exp2{t* -tD) = o(^-^^)^0{4>S), (28) 

where 4> — + (f)^)/2^, and S = — 0<|/2^. 
Within the volume V we do not discriminate between 
field amplitudes which differ by 0{ijl), and therefore take 
5 = 0(1). For (j) we note that, if to were to equal t* , then 
(jP' = 0{\/X) ^ 1, and in general (j) > 1. As a result, if 
there are no large numerical factors, we have 

tD<t*, (29) 

and the density matrix has become diagonal before the 
transition is complete. Detailed calculation shows jfl 
that there are no large factors [l^ . 

We already see a significant difference between the be- 
haviour for the case of a biquadratic interaction with an 
environment given in H1U|I and the more familiar linear 
interaction usually adopted in quantum mechanics. This 
latter would have replaced A//i^ just by incapable of 
inducing decoherence before the transition is complete. 

We note that, once the interaction strength is suffi- 
ciently weak for classical behaviour to appear before the 
transition is complete, this persists, however weak the 
coupling becomes. It remains the case that, the weaker 
the coupling, the longer it takes for the environment to 
decohere the system but, at the same time, the longer it 
takes for the transition to be completed, and the order- 
ing (|29|l remains the same. This is equally true for more 
general quenches provided the system remains approxi- 
mately Gaussian until the transition is complete. 



VI. EXTENSIONS OF THE MODEL 

Finally, it has to be said that taking only the short 
wavelength modes of the field as a one-loop system envi- 
ronment is not a robust approximation. We should sum 
over hard thermal loops in the 0-propagators. To be in 
proper control of the diffusion we need an environment 
that interacts with the system, without the system hav- 
ing a strong impact on the environment. We are helped 
in that, in the early universe, the order parameter field 
will interact with any field x for which there is no selec- 
tion rule. Again, it is the biquadratic interactions that 
are the most important. 

The most simple additional environment is one of a 
large number iV 3> 1 of weakly coupled scalar fields Xa, 
for which the action is extended to 



=^M+5[x] + 5int[0,x], (30) 



where S\(t)\ is as before, and 

N 

5i„t[0,x] = - E f / d^^4>\x)xl{x). (31) 

a=l 

where > 0. For simplicity we take weak couplings 
A ~ 5a and comparable masses ma — /i. The effect of a 
large number of weakly interacting environmental fields is 
twofold. Firstly, the Xa fields reduce the critical temper- 
ature Tc and, in order that — ^ ^ fi^, we must 

take A -I- X^ffa ^ 1- Secondly, the single x-loop contri- 
bution to the diffusion coefficient is the dominant x-field 
effect if, for order of magnitude estimates, we take identi- 
cal 5a = g/VN, whereby 1 > 1/VN » g ~ A. With this 
choice the effect of the (/)-field on the Xa thermal masses 
is, relatively, 0{1/\/N) and can be ignored. We stress 
that this is not a Hartrec or large-N approximation of 
the typ e that, to date, has been the main way to proceed 
[TsL fig . ITtI I for a closed system. 

Provided the change in temperature is not too slow the 
exponential instabilities of the (/)-field grow so fast that 
the field has populated the degenerate vacua well before 
the temperature has dropped to zero. Since the tempera- 
ture Tc has no particular significance for the environment 
field, for these early times we can keep the temperature 
of the environment fixed at = 0(Tc) (our calculations 
are only at the level of orders of magnitude). As before, 
we split the field as (/) = -|-<^> . The x-fields give an ad- 
ditional one-loop contribution to D{ko,t) with the same 
u{s) but a G++ constructed from (all the modes of) the 
X-field. The separation of the diffusion coefficient due to 
X into fast and slow factors proceeds as before to give 
a term that is identical to Hlt)|) (or (|23|l 'l but for its 
prefactor. 

Diffusion effects are additive at the one-loop level, and 
the final effect is to replace A^ in by + g"^ > A^, 
while leaving (|27|l unchanged. Although the relationship 
between Tc and A has been uncoupled by the presence of 
the Xa, the relationship l|28|) persists, with an enhanced 
right hand side, requiring that (|29|l is even better satis- 
fied. 

Given that the effect of further environmental fields is 
to increase the diffusion coefficient and speed up the onset 
of classical behaviour, additional fields interacting with 
the 4> field seem superfluous. However, the symmetries of 
the universe seem to be local (gauge symmetries), rather 
than global, and we should take gauge fields into account. 
We conclude with some observations from our work in 
progress |l8l | with local symmetry breaking. 

Local symmetry breaking is not possible for our real (j) 
field but, as a first step 0, it is not difficult to extend 
our model to that of a complex 0-field. At the level of 
0{2) global interactions with external fields and with its 
own short-wavelength modes, things are largely as before. 
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Local !7(1) symmetry breaking is most easily accommo- 
dated by taking the 0-iield to interact with other charged 
fields X through the local U{1) action 

S[cl>,A^,X\=S[<P,A^]+S^[A^,x]. (32) 
in which S[(j), A^] — 

(33) 

and 

S[A,„x]^ j d^x{{D^xrD^^X + ^\*x]- (34) 

For simplicity we have taken a single x held. The the- 
ory H32|l shows a phase transition at temperature Tc, and 
we assume couplings are such as to make this transition 
continuous. At the level of one loop the additional term 
to the diffusion function has derivative couplings. Hav- 
ing made a gauge choice, these give rise to explicit mo- 
menta factors in the generalisation of F. Unlike the 
contributions to D that we have seen so far, which are 
largely insensitive to the momentum scale fco, these con- 
tributions are strongly damped at large wavelength. In 
consequence, they barely enhance the onset of classical 
behaviour but, given that the effect of the other environ- 
mental modes is to enforce classical behaviour so quickly, 
it hardly matters. 



VII. CONCLUSION 

The previous paragraph says it all. For fast quenches 
weakly coupled environments make a scalar order param- 
eter field decohere before the transition is complete, un- 
der very general assumptions. An essential ingredient 
for rapid decoherence is nonlinear coupling to the en- 
vironment, inevitable when that environment contains 
the short-wavelength modes of the order parameter field. 
Had we onl y co nsidered linear coupling to the environ- 
ment, as in 19], for example (but an assumption that is 
ubiquitous in quantum mechanical models, from Brow- 
nian motion onwards) decoherence would not have hap- 
pened before the transition was complete, and we would 
not know how to proceed, although classical correlations 
would have occurred. 
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